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Bistatic  Scattering  Cross  Sections 
of  a  Composite  Rough  Surface 


1.  INTBODIICTION 

Perhaps  the  most  commonly  used  method  In  the  study  of  scattering  from 
randomly  rough  surfaces  Is  the  Rlrchhoff  or  physical  optics  method  [Beckmann 
and  Splzzlchlno,  1963;  Bass  and  Fuks,  1979]^*^.  Inherent  In  this  procedure  Is 
the  notion  of  Illuminated  region  and  shadowed  region.  However  It  Is  assuaied 
for  the  sake  of  simplicity  that  the  entire  surface  Is  Illuminated  [Stogryn, 
1967;  Barrlck,  1970]^’*.  Obviously  this  as8uaq>tlon  Is  Inappropriate  when  the 
angles  of  Incidence  and  observation  are  large.  Sancer  [  1969]^  has  considered 
this  problem  and  has  derived  *  shadow-corrected*  scattering  cross  sections. 
Several  others  [Brown,  1978;  Bahar,  1981  have  used  his  method  to  accoxint 
for  shadowing.  However  we  observe  that  his  method  leads  to  certain 
discontinuities  In  the  blstatlc  scattering  cross  sections .  This  Is  clearly 
unphyslcal.  Therefore,  In  this  report  an  expression  Is  derived  for  ‘shadow- 
corrected*  scattering  cross  section  that  does  not  suffer  from  the  above- 
mentioned  defect. 

The  contents  of  this  report  are  organized  as  follows .  The  geometry  of  the 
problem  Is  described  In  Section  2  vdille  the  statement  of  the  problem  Is  given 
In  Section  3.  In  Section  4  the  visibility  function  Is  derived.  Section  5 
contains  expressions  for  ‘shadow-corrected*  blstatlc  scattering  cross- 
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sections .  In  Section  6  we  have  a  brief  discussion  of  some  characteristics  of 


our  results .  The  conclusions  are  stated  In  Section  7 . 

2.  cnoMRur  or  not  problbi 

The  rough  surface  z  ■  ^  (x«y)  Is  represented  here  as  a  normally 
distributed  stationary  random  process  with  zero  mean  and  root  mean  sqtiare 
(rms)  height  .  The  mean  surface  coincides  with  the  xy-plane.  Let  T  and  R 
denote  the  locations  of  the  transmitter  and  receiver  and  let  0^  and  6^  be 
respectively  the  elevation  angles  of  IT  and  IR  (see  Figure  1) .  An  alternative 
coordinate  system  (u,v)  Is  chosen  for  the  xy-plane  such  that  the  projection 
(onto  the  xy-plane)  of  IT  coincides  with  the  u-axls  while  the  projection  of  IR 
coincides  with  the  v-axls .  The  angle  between  the  u-  and  v-  axes  Is  .  The 
slopes  of  IT  and  IR  are  denoted  as  p  and  v;  that  Is  to  say  ^  »  cot  6^.  >  v  »  cot  0^  . 
Let  p  m  and  q  m  denote  the  slopes  of  the  surface  at  an  arbitrary  point 

OU  O  V 

^(u«v) .  The  corresponding  rms  slopes  are  denoted  as  and  The  surface  Is 
assumed  to  be  Isotropic,  which  Inqilles  that  « ^2^ .  Further  for  a  normally 

distributed  surface  we  have  =  V2  2  'where  t  Is  the  correlation  length  of 
the  surface. 

3.  STATBBHT  OF  TBB  FROBUM 

It  Is  clear  that  when  the  angles  of  Incidence  and  observation  are  large 
the  phenomenon  of  shadowing  must  be  Incorporated  In  the  Klrchhoff  method. 
Sancer*s  shadowing  function  Is  Intended  to  serve  precisely  this  purpose.  Our 
objective  Is  to  obtain  'shadow-corrected*  scattering  cross-sections 
appropriate  for  blstatlc  situations . 
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The  problem  Is  to  find  the  probability  V(p,v|  p^,  q^)  that  an  arbitrary 
point  on  the  surface  Is  visible  to  both  transmitter  and  receiver. 

V  Is  widely  known  as  the  ‘shadowing  fimctlon’ .  But  It  Is  apparent  from  Its 
definition  that  the  above  label  Is  a  misnomer.  A  more  appropriate  name  for  V  Is 
'visibility  function* .  We  next  take  up  the  task  of  deriving  V. 

4.  VISIBILITY  yUWCTIOW 

It  Is  clear  that  for  the  point  IC^.p^.q^)  to  be  visible  to  both 
transmitter  and  receiver  the  preliminary  conditions  M'  >  p^  and  v  >  q^,  must  be 
satisfied.  Thus  If  the  surface  Intercepts  the  rays  IT  and/or  IR  It  must  do  so 
with  slopes  greater  than  that  of  the  rays. 

Consider  the  section  of  u>  and  v-axes  enclosed  In  the  circle  of  radius  a 
about  the  origin.  The  line  AOB  Is  divided  Into  2n  subsections  such  that  the 
odd-ntunbered  subsections  are  along  the  u-axls  while  the  even-nunibered 
subsections  are  along  the  v-axls  (see  Figure  lb) . 

We  denote  the  event  that  the  Interception  of  one  of  the  rays  occurs  within 
the  subsection  1  as  The  cosqplement  of  the  event,  namely,  the  event  that 
neither  of  the  rays  Is  Intercepted  within  the  subsection  1,  Is  denoted  as 
Consider  ■  P(H  .  .  .  ^52^) .  This  denotes  the  probability  that  there 

Is  no  Interception  of  the  rays  by  the  surface  within  the  radius  a.  can  be 

expanded  as  follows . 
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V 

a 


1  - 


I 


i-1.3 


P(H^) 


I 

1-2.4 


P(Hp 


+  "  0.5  P(H.n  H.) 

i.j-1.3  ..  ^  J 


+  "  0.5  P(H.n  H.) 

i.j-2.4  ..  ^  J 

i'-j 


I 

j-2.4 


P(Hj^n  Hj)  -  . 


In  the  limit  as  n  ->  <»  and  as  the  subsections  become  uniformly  small 
we  have 


(1) 


V 

a 


1  -  J  du.  g(u. 
0 


a  a  a 


a  a  a  a 

+  .5  /  dv.  /  dv«  g(v.,,v,)  +  /  du.  /  dv.  g(u,  ,v  )  -  .  .  . 
0^0^  ^  ^  0^0  ^ 

(2) 

where  the  conditional  density  functions  g*s  are  derived  In  Appendix  A. 

Our  visibility  function  V  Is  now  readily  obtained  as 

V  -  lim  V 

a-»  •  ®  (3) 


But  Instead  of  using  Eq.  (2)  we  use  an  alternate  representation  that  Is 
convenient  for  our  later  approximations . 
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«0  <0  «0  CO 

+  .5  J  dv.  /  dv-  h(v-,v-)  +  /  du  /  dv,  h(u.  ,v  ) 
0'*'0^  0^0"^ 


where  h’s  are  the  ctumilant  functions  [Rubo,  1962] All  the  cumulant  functions 
can  be  expanded  In  terms  of  density  functions.  For  exa]iq>le, 


-  g(Uj^) 

-  g(v-) 

^  ^  (5) 

h(u^,U2)  -  g(uj^,U2)  -  g(Uj^)  gCUj) 

We  observe  that  when  shadowing  Is  significant  the  distances  between  crossing 
points  become  large.  That  Is  to  say,  |  Uj^  -  uj  |  >  In  otherwords  adjacent 
crossing  points  become  uncorrelated.  Further,  an  l]q>ortant  property  of  the 
higher-order  cuanilant  functions  Is  that  they  vanish  rapidly  as  the  associated 
random  quantities  tend  to  become  uncorrelated.  Thus  to  a  first-order,  Eq.  (4) 
may  be  approximated  as 


V  -  exp 


{•  «o  «o  CO  ] 


(6) 


FromEqs.  (A3)  and  (5)  we  have 


h(uj^)  -  /  dp  (p-p)  f({^+/iuj^,p) 

M 

We  further  note  that  the  slopes  and  height  at  any  point  on  the  surface  are 


(7) 
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uncorrelated.  This  iii?>lies  that  f  (4Q+fiuj,p)  =  f  (^q+huj)  f  (p) .  Thus  on 
evaluating  Eq .  ( 7 )  we  obtain 

h(u. )  -  f «  +/XU  )  A.  (n) 

1  oil  (8a) 


i^ere  ^ 

-  77?^  ■  "f-  [  jfl]  '**'> 

Similarly  we  obtain  h(v^)  as 


h(V^)  -  f(C^+  i'^y)  A^(i/)  (9) 

We  note  fromEq.  (5)  that 

h(Ui,Vi)  -  g(Uj^,Vj^)  -  g(u^)  g(Vj^)  (10) 


vdiere  gCu^.v^)  is  defined  in  Appendix  A  as 
00  00 

-  /  dp  (p-m)  J  dq  (q-i/)  (11) 

ft  V 

In  principle  f  Eq .  (11)  can  be  evaluated  using  the  general  expression  for  f  ( X ) 
as  derived  in  Appendix  B.  But  it  appears  too  con^licated  to  obtain  analytical 
results .  We  therefore  focus  attention  on  two  special  cases . 


Case  1  :  Large  4>:  >  1 

In  this  case  r/^j  »  1,  and  as  shown  in  Appendix  C,  f(X)  assumes  the 
simplified  form 


f(X)  -  f(Xj^)  f(X2)  f(X3)  f(X^) 


(12) 
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Hence 


g(u^.Vi)  -  /  dp  (p.#t)  f  dq  (q-i/) 

M 

.  f(p)  f(«^+»/v^)  f(q) 

FromEqs.  (13) ,  (10),  (8),  and  (9) 


h(Ul,Vi)  -  0 


Thus  the  visibility  function  in  this  case  is 


Using  Eq.  (8)  we  have 


/  du^  h(\ij^)  -  J 


”2jr“ 


erfc 


Likewise 


/  dv  h(v  ) 
0  ^  ^ 


r  f  1 

IT"  -7^ 


(13) 


(14) 


(15) 


(16) 


(17) 
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Thus 


exp 


_L 

2 


erfc 


where 


(18) 


(19a) 


Aj^(*') 


Aj^(i/) 


(19b) 


As  explained  earlier,  V  represents  the  conditional  probability  that  a  point 
on  the  surface  with  height  ^  and  slopes  p^,  is  visible  to  both  transmitter 

and  receiver.  We  also  mentioned  that  to  avoid  self- shadowing  the  inequalities 
pQ  <  p  and  q^  <  V  must  hold .  The  above  ideas  are  condensed  in  the  following 
expression  for  V,  or  more  specifically  V(p.,vl  ^^jp^^q^) . 


M-Pq) 


•  exp 


erfc 


7?^J  (  ^i^^>  ^i<*'>/ 


(20) 

where  X  (t)  is  the  Heaviside  function  with  the  usual  definition  that  X  (t)  =  0 
for  t  <  0  and  X  ( t )  =  1  for  t  S  0 .  From  Eq .  (20)  we  can  obtain  the  visibility 
function  V(n,v I  Po^q©^ 
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(21) 


In  this  case  f  (X )  takes  the  form  (see  Appendix  C) 


f(X)  -  f(X^)  f(X2)  SCXj^-Xj)  fi(X2-X^) 

so  that 

6(u^’V^)  f(p-q) 


(25) 


Hence 


g(Ui,Vi)  -  /  dp  (p-m)  f  dq  (q-i/) 
M 


•  f(p)  «(p-q) 


f(^^+MU^)  5(u^-v^)  f  dp  (p-M)^  f(p) 

M 


where 


A2<At) 


«p  [ - -^ ]  { - 1  * y('/2)  4  [  jfi] } 


+  0.5  erfc  [  1 

72 


It  follows  fromEqs.  (10),  (26),  (16),  and  (17)  that 


/  du.  /  dv  h(u.  ,v.) 
0^0  ^ 


«0  OO 


-  J  -  g('ii)  gCvj^)  ] 
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to  <D 


/  du^  J  dv^  ^2^^^ 


-  /  dxi^  g(Uj^)  J  dv^  g(v^) 


-|-  A^itt)  erfc 


75^J 


1—9  2 

A^(m)  erfc 


J7A 


IJ 


(28) 


where 

A2<m)  -  — —  (29) 

FromEqs.  (28),  (16),  (17) ,  and  (6)  we  have 


/  -2 

V(/i,MKo.Po)  -  exp  -  Aj^(m)  erfc 


-  [  2  Aj^(m)  -  A2(m)  ]  erfc 


(30) 


The  visibility  function  independent  of  is  given  as 
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V(M.M|p„)  -  /  f(«o) 


*o  •  o 


-  ^ 


4a 


exp 


4a 


erf 


2a 


-  erf 


(31) 


where 


a  -  Aj^(m) 


(32a) 


-  2Aj^(m)  -  AjCm) 


(32b) 


As  explained  In  case  1  the  above  expression  assumes  that  p  >  Pq  and  hence  It  Is 
Identical  to  V(p,p|Pq). 


5.  scATmtzno  citoss  sscrzom 


Analytic  solutions  for  scattering  cross  sections  are  available  at 
present  only  for  certain  types  of  random  surfaces.  The  one  among  them  of 
particular  Interest  to  us  here  Is  the  composite  rough  surface  [Beckmann, 
1965;  Bahar,  1981]^’^  .  Both  for  Its  generality  and  Its  utility  In  various 
applications,  the  composite  rough  surface  deserves  attention.  As  the  name 
suggests ,  this  Is  composed  of  a  ssiall  scale  rough  surface  superposed  on  a 
large  scale  rough  surface.  We  characterize  the  large  scale  roughness  by  rms 
height  bi  and  correlation  length  .  The  corresponding  parameters  for  small 
scale  roughness  are  ^2  and  /2*  practical  applications  where  the  composite 
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rough  surface  model  Is  useful  It  Is  found  that  the  two  scales  are 
uncorrelated.  This  Is  primarily  because  the  large  scale  roughness  and  the 
small  scale  roughness  are  usually  generated  by  different  physical 
mechanisms.  In  this  case  It  has  been  shown  [Semenov,  1966;  Barrlck  and  Peake, 
1968  ]  that  the  problem  can  be  broken  up  Into  two  parts .  Scattering  due  to 

large  scale  roughness  Is  determined  using  the  Klrchhof  f  method  [Beckmann  and 
Splzzlchlno,  1963]^  while  scattering  due  to  small  scale  roughness  Is 
evaluated  using  the  perturbation  method  [Rice,  1951]^^.  The  above  two 
quantities  are  then  added  Incoherently  to  obtain  the  scattering  due  to  the 
composite  rough  surface . 

It  should  be  mentioned  that  the  above  procedure  Is  valid  only  as  a  single 
scattering  approximation.  Thus  we  agree  to  Ignore  multiple  scattering 
effects .  Still  there  Is  another  factor  to  be  taken  into  consideration  -  the 
phenomenon  of  shadowing.  Sancer  [  1969]^  has  shown  that  for  large  scale  rough 
surfaces  the  *  shadow- correction*  Is  incorporated  by  multiplying  the  usual 
scattering  cross  sections  by  the  visibility  function  V  (m.,v)  .  Further,  Brown 
[  1978]^  has  formally  derived  the  scattering  cross  sections  of  a  composite 
rough  surface  by  a  perturbation  scheme  to  deal  with  the  small-scale 
fluctuations  of  the  underlying  large  scale  rough  surface.  This  leads  to  the 
conclusion  that  the  *  shadow-correction*  to  the  composite  rough  surface  may  be 
obtained  similarly  by  using  the  saaie  visibility  function.  In  other  words,  the 
shadow- corrected  scattering  cross  section  of  the  composite  rough  surface  Is 
written  as 


T^cd“ 


\  c.d  /  -  \  V,H  / 


(33) 


where  V  (^tV)  Is  the  visibility  fimction  derived  in  Section  4.  The  scattering 
cross  sections,  and  ate  due,  respectively,  to  the  small  and  large 

scale  rough  surfaces  that  constitute  the  given  coiiq>oslte  rough  surface. 
Barrlck  and  Peake  [  1967]  “  have  calculated  and  .  .  Their  results  are 

quoted  here . 


where 


,2  r  1  .2  ,2  ,  2  ^  2.  1 

-  *  *2  •*>’1  ■  »  “o  ‘2  <’x  ’y>  J 


ij  -  sin  t.  -  sin  $  cos  4 
X  i  s 


(36a) 


»  -  sin  $  sin  4 

y  s 


(36b) 


[  +  (/*y-l)(MjSintfj^sinJ^-  tj^t^cos^)] 


(37a) 


(t.+  «  cos#,)  (t  +  u  cos#  ) 
1  r  1  s  r  s 


ti.r(Mr-l)  - 
(tj^+P^COS#  J^)  (t^+c^cos#^) 


(37b) 


t  n  (c  -1)  —  t.c  (fi  -1) 
s^r'  T  '  i  r''^r  '  .  , 

-  sinf 

(t^+/i^COS#^)  (tj^+C^C0S#j^) 


(37c) 
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(37d) 


\v 


[  *r(My-l)cos^  +  (c^-l)(c^sin0^siTi$^-  tjt^cos^)J 


(t.+  £  costf.)  (t  -»•  £  costf  ) 
1  i  X  s  r  s 


r  2  1 

-  [  Vr  -  ^  ] 


.  X  -  i.s 


(ii) 


'cd 


-  I  h/ 


where 


r(*) 


*2  2 


exp 


J  \  . 


2  .  2 


\ 

/ 


1?^  -  -  COsSj^-  cos^ 


B, 


HH 


"  "  *2*3  ^  ^ 


®VH  ”  "4^  [  ■  ^^  *3  ^  ^  *2  ^  ] 


B 


HV 


“  *2  ^  -  *V  *3  ^  ] 


V  ■  ■  [  *Sl  ^  ^  +  Rv  *2*3  ] 


(38) 


(39) 


(40) 


(41) 


(42a) 


(42b) 


(42c) 


(42d) 
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where  Rg  and  Ry  are  the  Fresnel  reflection  coefficients  with  angle  of 
incidence  y  (defined  below)  and 


1  +  sin  $.  sin  6 
1  s 


cos  ^  -  cos  cos  $ 


(A3a) 


Sj  “  cos  sin  0^  +  sin  8^  cos  8^  cos  ^ 


(A3b) 


a,  -  sin  8.  cos  8  +  cos  8.  sin  8  cos  d 

3  IS  1  s 


a,  »  cos  8.  +  cos 
4  i 


(43c) 


(43d) 


^  «■  arccos\.5  [  1  *  sin  8.  cos  8  cos  ^  +  cos  8.  cos  6 
L  X  S  X  2 


]) 


.5 


(44) 


The  iiq>ortant  symbols  in  the  above  expressions  are  defined  next.  The 
elevation  angles  of  the  incident  and  the  scattered  rays  are  9^^  and  9^.  They  are 
identical  to  6^  and  9^.  mentioned  in  the  previous  sections  of  this  report .  The 
propagation  constant  in  free  space,  the  condition  above  the  rough  surface  is 
kg.  The  mediiim below  the  rough  surface  is  of  relative  permittivity  and 

relative  permeability  The  Fresnel  reflection  coefficients  for  horizontal 
polarization  Rg  and  for  vertical  polarization  Ry  are  given  as 


2  5 

cosy  -  -  sin  y  )' 

2  5 

cosy  +  -  sin  y  )■ 


(45a) 
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(45b) 


2  5 

cos^  -  -  sin  )■ 

“  2  T 

cos^  +  (c^j.  -  sin  0  )* 


6.  DXSCnSSKW 


Our  interest  here  is  to  study  the  expressions  for  the  shadow-corrected 
bistatic  scattering  cross-sections  given  in  Eq.  (36) .  Barrick  and  Peake 
1 1967]  have  elaborately  analyzed  the  characteristics  of  and  ;  so  it 
suffices  for  us  here  to  study  the  behavior  of  the  visibility  function  V  (ji,v) . 
We  notice  that  V  (p,v)  has  a  rather  complicated  functional  dependence  on  the 
physical  parameters  p,  v,  6^,  6^  and  <|i.  therefore,  we  must  be  content  with 
looking  at  two  limiting  cases . 

As  the  visibility  parameters  A  *  ^  ® 

clear  that  the  surface  area  of  visibility  diminishes,  'uhlcb  implies  that  V  (p, 
v)  shotU.d  vanish  in  the  above  limit.  In  order  to  verify  this  physical  fact  we 
proceed  as  follows.  On  using  the  asynqitotic  expression  erfc(t)  -  I  -  Ztl-Jn 
for  t  «  1  we  obtain  fromEq.  (23)  the  result. 


V(n,v)  -  Zjn  — 

p  +  1/ 


Ji  ,  P  <1 


We  thus  see  that  V  (p,v)  linearly  approaches  zero  as  jl,  v-^0  . 

On  the  other  hand  p  ,  v  »  1  when  we  use  the  asymptotic  expression  erfc(t)  - 
exp(-t^)  /Vit  for  t  »  1  and  find  that  V  (p,v) approaches  unity.  This  is  an 
intuitively  satisfying  physical  result.  For  in  this  limit  virtually  the 
entire  surface  is  visible  and  hence  there  is  no  need  for  *  shadow-  correction* . 
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It  Is  only  In  this  region  that  the  restilts  of  Barrlck  and  Peake  [  1967]  are 
strictly  valid.  Our  results  are  hence  the  natural  generalizations  of  theirs . 

It  Is  apparent  that  our  shadow- corrected  blstatlc  scattering  cross- 
sections  have  several  other  Interesting  characteristics.  But  since  their 
analytic  forms  are  very  complicated  It  seems  that  they  can  be  further  studied 
only  by  considering  some  numerical  examples. 

7.  COMCLUSION 

We  have  considered  a  normally  distributed  composite  rough  surface  and 
derived  the  shadow-corrected  scattering  cross  sections .  This  contrasts  to 
several  other  works  where  the  notion  of  shadowing  Is  Ignored.  A  key  role  Is 
played  In  the  analysis  by  the  visibility  function.  Physically  this  function 
represents  the  fraction  of  the  surface  visible  to  both  the  transmitter  and  the 
receiver.  A  detailed  derivation  of  the  visibility  function  Is  provided.  Also 
Included  Is  a  discussion  of  the  properties  of  the  visibility  function.  This 
work  stay  be  regarded  as  an  la^rovesient  of  the  work  of  Sancer  [  1969]^  whose 
expressions  for  scattering  cross  sections  contain  certain  erroneous 
discontinuities .  This  work  also  has  taportant  practical  applications . 
Several  radars  presently  must  operate  at  grazing  angles .  Shadowing  Is  a 
dominating  physical  phenosienon  In  such  situations .  Since  many  modem  radars 
operate  In  the  blstatlc  mode  and  since  the  conq>oslte  rough  surface  Is  a  very 
useful  practical  siodel  we  anticipate  that  our  studies  In  this  report  will 
serve  to  enhance  the  performance  capabilities  of  such  radars . 


Appendix  A 

Conditional  Density  Functions 


According  to  Its  definition,  g(u^)Auj  represents  the  conditional 
probability  that  the  surface  intercepts  the  ray  IT  in  the  interval 
(u^,  uj-t-Auj)  with  p>n  given  that  the  height  at  I  is  and  its  slope  is  p^. 
Suppose  the  surface  crosses  IT  at  u  in  (Uj^,  uj^+Au^) .  Then  at  the  crossing 
point,  §Q+fiu  =  §(ui)+p(u-Uj) .  This  inplles  that  >  ^(ui)  >  - 

(p-p)AU|  .  Let  the  event  that  the  surface  crosses  IT  with  p  >  p  in  the  interval 
(uj^,  Ujj^'t'AU]^)  bed.  Then 


•  ^  +/iU, 

P(^  )  -  J  dp  J  **  ^  f(?.P) 


-  AUj^  J  dp  (p-m)  f(?Q+MUj^.P) 

where  the  approximate  result  is  obtained  using  the  mean  value  theorem. 
Throughout  this  report  we  denote  normal  density  functions  by  f '  s .  For  exaiiq>le 
f  (^tP)  denotes  the  joint  normal  density  function  of  slope  and  height  at  u. 

Now  if  the  point  I  from  which  the  ray  emanates  is  of  height  ^  and  slope  p^ 
we  can  deduce  the  conditional  probability  g2^(uj)Au]^  fromEq.  (Al)  as  follows. 


g(uj)Au^  -  fTiTTT  ^ 

^0*^0  /i 
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When  shadowing  Is  significant  it  Is  apparent  that  the  distances  between 
adjacent  crossing  points  are  large.  This  means  that  adjacent  crossing  points 
are  weakly  correlated.  Then  we  can  approximate  f(^+pu2,  p:^q,Pq)  as 
fC^+pupp)  fC^.pp).  Thus  we  have  fromEq.  (A2) 


w 

/  dp  (p-M)  p) 


Similarly  we  can  proceed  to  obtain  other  conditional  density  functions .  For 
example, 


g(Vj^)  -  /  dq  (q-i/)  f(^^+i/Vj^,  q) 


g(Uj^,Vj^)  -  /  dp  (p-p)  f  dq  (q-»/)  f (^^+AiUj^,p;^^+>/v^,q) 
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Appendix  B 

The  Density  Function  f  (^(u,0),p;^(0,v),q) 


Let 


-  ^(u.O) 


(Bla) 


*2  -  IT  ” 


(Bib) 


X3  -  ?(0.v) 


(Blc) 


X4  -  iT  "  ** 


(Bid) 


How  our  object  in  this  appendix  is  to  cosqmte  the  fourth-order  normal  density 
function  f  (Xi,X2,X3tX4)  or  f(X) .  In  terms  of  its  characteristic  function  we 
have 


f(X) 


*(W)  exp[-  iX*W  ] 


(B2) 


Hoting  that  X^ ,  X2 1  X3 ,  X4  have  zero  means  we  have 


#(ff)  -  exp  f“  iW«K»W  ] 
where  K  is  the  covariance  matrix  defined  as 


(B3) 


-  cov  (Xj^,  Xj)  “ 


(B4) 
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Since  the  surface  is  assumed  to  be  statistically  homogeneous  and  isotropic 
the  height-height  correlation  function  C(r  )is  given  as 

C(r)  -  CCIrj^-rjD-  (  (B5) 

In  terms  of  C(r  }  the  covariance  matrix  K  takes  the  form. 


C(0) 

0 

C(r) 

C  (r) 

V 

0 

-  C"(0) 

C  (r) 
u 

C  (r) 
uv' 

C(r) 

C  (r) 
u 

C(0) 

0 

C  (r) 

C  (r) 
uv'  ■' 

0 

-  C"(0) 

(B6) 

where 


C.(,) 

dr 


c  (r)  . 

u  au 


C  (r)  .  ^ 
V  av 


C  Cr)  - 
uv^  ^  au  av 


(B7> 


2  2 

-  u  +v  —  2uv  cos^ 


(B8) 


For  definiteness  and  for  the  purpose  of  discussion  we  assume  that  the  surface 
heights  have  the  following  correlation  f^lnction 


2 


2 


C(r)  -  exp  (  -  ] 


In  this  case  the  various  elements  of  the  covariance  matrix  are  given  as 
follows . 


C(0)  - 


C-(0)  -  -  < 


C^(r)  ~  (U'V  cos^)  exp  (  -  '' 

C^(r)  -  -  (v-u  cos^)  exp  (  -  ) 


(BIO) 


(BID 


(B12) 


(B13) 


C  (r)  - 
uv' 


-  aj  exp  (  -  (r/t^)^  ) 


/  2  \ 
\  cos4  +  — ^  (u-v  cos^)  (v-u  cos4)  f 


(B14) 
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Appendix  C 

Density  Function  For  Two  Special  Cases 


Case  1  :  large  i,  i/ty>  1 

FromEqs.  (B9)>(B14)  we  note  that 


C(r)  -  C^(r)  -  C^(r)  -  C^(r)  -  0 


Thus 


K 


The  characteristic  ftmctlon  Is  then 


(Cl) 


(C2) 


♦(W)  -  exp  {  -  }  (C3) 

Substituting  Eq.  (C3)  IntoEq.  (B2)  we  have 

f  (X) - ^  f  d^W  exp  {  -  (\\^  +  ^1^3^  ^ 

(2x)  — 

.  exp  [  -  1  X2W2+  X3W3+  X^W^)] 

(CA) 
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where 


40 

h  -  IT  J  <Wi  exp  [  -  .5  -  i  XjWj  ] 


40 

h--^  ;  dWj  .xp  [  - 


h--^  fdWj  exp  [  - 


i  XjWj  ] 


40 

^4  -  2^  /  dw,  exp  [  -  .5  /.^  -  1  ] 


These  integrals  are  elementary  and  the  results  are 


^1  “  7T2^ 


1  r  ^1  1 

[  -  ^  J  -  ^(^i> 


f(X2) 


f(X3) 


fCX^) 


(C9a) 


(C9b) 


(C9c) 


(C9d) 
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(CIO) 


f(X)  -  f(X^)  fCX^)  f(X3)  f(X^) 


Case  2:  i"0:r-0:u»i/ 


In  this  case 


C(0)  -  a! 


C  (0)  -  C  (0)  -  0 
u  V 


C  (0)  -  < 

UV  ■'  1 


Thus  the  covariance  matrix  takes  the  following  form 


(Cll) 


The  four  eigenvalues  of  K  are  2b^  ,  21-^  »  0,  0  .  Corresponding  to  them  we  may 
choose  four  linearly  independent  eigenvectors  as 
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Let 


(C12) 


Note  that  we  have  normalized  the  eigenvectors .  By  a  well  known  theorem 
in  matrix  analysis , 


K  f  -  A 


where 


diag  (  2a2  ,  2^2  ,  0  .  0 


Also  we  know  that  normalized  eigenvectors  form  an  orthonormal  set . 
This  implies 


f  f '  -  f  - 1 


But  since  in  our  case  P  -  P  ^  we  have 


P  -  P 


The  density  function  f  ( X )  is  given  as 


f(X)  -  (2jr)'^  J  d^W  exp  [  -  0.5  W.K»W  -  iX*W  ] 


On  making  the  transformation 


W  -  p.z 


(C13) 


(CIA) 


(C15) 


(C16) 


(C17) 


(CIS) 
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(C19) 


A— 

d  W  -  J  d  Z 


where  J  is  the  Jacobian  of  the  transformation.  In  our  case 


J  -  det  P  -  1 


f(X)  -  (2jr)*^  /  d^Z  exp  [  -  0.5  Z»P  K  P»Z  -  iX»P*Z  ] 


(C20) 


-  (2ir)  ^  f  d^Z  exp  [  ~  0- 


5  Z«A‘Z  -  iX-P-Z 


(C21) 


0.5  Z«A«Z 


-  a2  z2  +  -I  ^2 


(C22) 


X»P«Z  - 


^  {  X^(Z^+Z3)  +  X2(Z2+Z^)+  X3(Z3-Z3)+  X^(Z2-Z^)  } 


(C23) 


Therefore 


f(X)  -  J2  J3  \ 


where 


(C24) 
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1  * 

■'i  ■  ST  /  “i  4  ] 


1  * 

■  ^  i  *^2  **P  [  *^1  4  ■  i(X2+X^)Z2/y2  ] 


1  * 

“  ^  /  <iZ3  exp  [  -  i(X^-X3)23/y2  ] 


1  * 

\-2r  /  <^4  •’T  [  -  ] 


All  the  above  integrals  can  be  readily  evaluated  and  the  results  are 


(C25a) 


(C25b) 


(C25c) 


(C25d) 


••i  ■  JuJ)  /.J  [  - 


■  ■  TT  ^<*1> 


■^2  -  7r«v 


J3  ~  J7 


~  J2  SCX^-X^) 


(C26) 


{C27) 


(C28) 


(C29) 


So  finally 


f(X)  -  f(X3^)  f(X2)  5(X^-X3)  «(X2-X^) 


(C30) 
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Nomenclature 


^<x.y) 


correlation  function 


normal  density  function 
conditional  density  function 


cumulant  function 


large  scale  ms  height  of  the  surface 
small  scale  rms  height  of  the  surface 


covariance  siatrlx 


large  scale  correlation  length  of  the  surface 
small  scale  correlation  length  of  the  surface 
slopes  along  the  transmitter  and  receiver  directions 
Fresnel  ref  lection  coefficient  for  vertical,  horizontal 
polarization 

blstatlc  scattering  cross  section 
large  scale  rms  slope  of  the  surface 
small  scale  rms  slope  of  the  surface 

slope  of  the  ray  emanating  from  I  towards  the  transmitter 
slope  of  the  ray  emanating  from  I  towards  the  receiver 
random  function  describing  the  z -coordinate  of  the  surface 
azimuthal  angle 
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